Dimensional Evolution of Spin Correlations in the Magnetic Pyrochlore Yb2Ti207 
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The pyrochlore material Yb2Ti207 displays unexpected quasi- two- dimensional (2D) magnetic 
correlations within a cubic lattice environment at low temperatures, before entering an exotic dis- 
ordered ground state below T=265mK. We report neutron scattering measurements of the thermal 
evolution of the 2D spin correlations in space and time. Short range three dimensional (3D) spin 
correlations develop below 400 mK, accompanied by a suppression in the quasi-elastic (QE) scat- 
tering below ^0.2 meV. These show a slowly fluctuating ground state with spins correlated over 
short distances within a kagome-triangular-kagome (KTK) stack along [111], which evolves to iso- 
lated kagome spin- stars at higher temperatures. Furthermore, low-temperature specific heat results 
indicate a sample dependence to the putative transition temperature that is bounded by 265mK, 
which we discuss in the context of recent mean field theoretical analysis. 

PACS numbers: 75.25.-j, 75.40.Gb, 75.40.-s, 75.10.Jm 



I. INTRODUCTION 

The corner-sharing tetrahedral geometry of the py- 
rochlore lattice favors strong geometric frustration^ 
wherein a competition of exchange interactions between 
localized magnetic moments results directly from their 
arrangement within the lattice. This phenomenon tends 
to suppress transitions to magnetic long range order 
(LRO), leaving in their place unusual short-range mag- 
netic correlations at low temperatures. The form of 
such correlations varies between frustrated systems, with 
ground states typically determined by subleading terms 
in the relevant Hamiltonian, such as dipole interactions, 
further neighbor exchange, and structural distortions 
within the underlying crystal lattice^ii^ The rare earth 
titanate series of materials, which have the chemical for- 
mula i?2Ti2 07 with R being a trivalent magnetic rare 
earth ion, crystallize into the face centered cubic (fee) 
space group Fd3m which characterizes the pyrochlore 
lattice. This series provides excellent examples of the va- 
riety of exotic ground states that can emerge in nearly 
identical chemical environments when geometric frustra- 
tion is at play. Several of the rare earth titanates exhibit 
disordered ground states which persist down to 30mK, 



namely the spin ices Dy2Ti207 and Ho2Ti207r'— and 
the enigmatic spin liquid Tb2Ti207^ii Others in the se- 
ries do order magnetically, such as Gd2Ti207^— and 
Er2Ti207r^ii^ but the nature of the ordering is unusual 
in both cases. In all of these cases, the source of the frus- 
tration (or lack thereof) is well-understood in terms of 
the spin interactions and the single-ion anisotropy im- 
posed by the crystal electric field. In particular, the spin 
ices combine effective ferromagnetic (FM) exchange with 
strong Ising-like anisotropy, such that the moments are 
constrained to point into or out of each tetrahedron. 

Yb2Ti2 07 provides an intriguing contrast to the spin 
ices. It is known to combine FM exchange (^ciy^SOOmK 
to SOOmKi^J^) with XY anisotropy J^-i^ The continu- 
ous degree of freedom provided by the XY anisotropy 
leads one to naively expect an ordered FM ground state. 
However, despite what appears to be a transition near 
200mK^i2Ai^i2£ most experimental evidence suggests that 
Yb2Ti2 07 displays a disordered ground state down to 
30mK in zero-magnetic field J^ii^ Application of a mod- 
est magnetic field along [110] at 30 mK induces a transi- 
tion to a polarized, ordered phase with well defined spin 
wave excitations.^^ Using anisotropic exchange Hamilto- 
nians, recent theoretical studies have sought to under- 
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FIG. 1: (color online) a) The temperature dependence of neu- 
tron scattering intensity at two Q points that lie on the rod 
of scattering, Q1.54 (red x's) and Q1.07 (blue o's) (lOK back- 
ground subtracted). Bottom Inset: the position of the two Q 
points is shown relative to the fee Brillouin zone boundaries 
in the (HHL) plane. Top Inset: the difference of intensities, 
Ii.or - I1.54, with the solid line as a guide to the eye. b) Elastic 
time-of-flight neutron scattering data at 30mK and 500mK 
(reproduced from Ref . \13i) . Error bars represent ± la. 



stand the spin wave dispersions in the magnetic field in- 
duced ordered stated as well as the diffuse scattering at 
relatively high temperatures in zero field^^ and the local 
spin susceptibility J^i^ Interestingly, some of these stud- 
ies indicate that the dominant exchange interaction is 
the Ising component, which may explain the lack of sim- 
ple FM behaviour in Yb2Ti207, making it an exchange 
analog to the spin icesJ^i^ 

An early heat capacity (Cp) measurement showed an 
anomaly at T=214mK, indicating a possible phase tran- 
sition, in addition to a broad Schottky-like hump at 
2K^ The putative transition, whose temperature varies 
somewhat in the literature and which we henceforth 
label Tc, has previously been explored through sev- 
eral techniques ^i^i^i^ Much of this characterization 
is consistent with the absence of conventional LRO be- 
low Tc, although single crystal neutron diffraction and 
AC susceptibility results presented in Ref. [25|, suggest 
a colhnear FM ground state. Hodges et al found a dis- 
continuous change in the spin fluctuation frequency at 
Tc=240mK as measured by Mossbauer and /iSR, along 
with an absence of magnetic Bragg peaks in neutron 
powder diffraction^ The specific heat of the sample 
used in Ref. |l^ was found to display a sharp anomaly 
at 250mK^ close to the observed first order drop in 
the spin fluctuation rate. Single crystal neutron scat- 
tering studies revealed a pattern of diffuse scattering, 
notably taking the form of rods of scattering along the 
<111> directions, with a characteristic QE energy scale 



of 0.3meVJ^i2S This diffuse scattering, which indicates 
short range spin correlations, is present both above and 
below Tc, but a qualitative difference was observed be- 
tween 500mK and 30mK that indicated short-range 3D 
spin correlations at the lowest temperatures J^ The ob- 
served change was that the rod of scattering became 
more strongly peaked near the (111) zone center, moving 
from a featureless rod at 500mK to a rod with structure 
at 30mK. The conventional interpretation of a flat rod 
of scattering in any diffraction pattern is the presence of 
2D correlations. Though the rod of scattering is flat at 
500mK, this interpretation is striking in the present case, 
as it implies a magnetic decomposition of the 3D py- 
rochlore structure into a stacking of alternating kagome 
and triangular planes along <111>, breaking the under- 
lying cubic symmetry. 

It should be noted that in Ref. '22', where exchange pa- 
rameters were extracted by fitting a modified measure- 
ment of the diffuse scattering pattern in the HHL plane 
at 1.4K , the real space spin correlations derived from 
these parameters were reported to show a more nuanced 
anisotropic behavior. However, their calculated spin cor- 
relation functions display some decoupling of the kagome 
planes with pronounced correlations along near neighbor 
[110] directions. 

The nature of the transition near 200mK, and the re- 
sulting ground state of Yb2Ti207, remains a matter of 
debate. Most of the experimental studies indicate a dis- 
ordered ground state. This is supported by a recent theo- 
retical study^^ that makes use of the measured high-field 
spin wave spectrum to determine microscopic exchange 
parameters. The parameters are highly anisotropic, fa- 
voring a local <111> exchange. From a mean field per- 
spective, these parameters predict an ordered ferromag- 
netic state below 3.2K. Frustration and quantum fluc- 
tuations therefore must be responsible for suppressing a 
transition to long range FM order, possibly to zero tem- 
perature. Furthermore, the extracted exchange parame- 
ters may place Yb2Ti207 within the non-magnetic spin 
liquid ground state studied by Refs. 27 and i28i. In this 
scenario, the zero-field QE diffuse scattering measured 
at 30mK in the most recent single crystal work, Ref. [l3|, 
would represent two-spinon-type magnetic excitations. 

Contrary to this scenario, the work of Yasui et al 
indicates an ordered, collinear FM ground state below 
250mK, albeit possibly a glassy one with a very long re- 
laxation time^ In light of the results to be described 
in the present report, which identify a prevalent sample 
dependence to the temperature and sharpness of the spe- 
cific heat anomaly at "Tc" , this FM ordering remains a 
possibility for the true ground state of Yb2Ti207. How- 
ever, we will show that the sample dependence of the 
Cp anomaly is seemingly irrelevant to the changes in 
energy-integrated diffuse scattering, which occur at a 
temperature that does not coincide with any specific heat 
anomaly. Furthermore, the spread of Tc in various sam- 



pies is quite small compared to the overall suppression 
from the predicted mean field transition temperature, 
lending further support to the idea of a low temperature 
state dominated by quantum fluctuations. 



a) 1600 



60mK 
300mK 
540mK 




-0.4 -0.2 0.2 0.4 
[-2AH,AH,AH] 



0.4 -0.2 0.2 

[-2AH,AH,AH] 



FIG. 2: a) (lOK background subtracted) Longitudinal rod 
scans at T = 60 mK and 300 mK. Data points at 540mK are 
consistent with a Q dependence that follows the Yb^^ mag- 
netic form factor (black solid line) . b) Examples of transverse 
scans across the rod, centered at Qi.o? and Q1.54. The insets 
show the positions of the scans in reciprocal space. 

This paper details the temperature dependence of the 
aforementioned changes in diffuse magnetic scattering in 
zero-fleld, and determines the extent to which they corre- 
late to features in the speciflc heat. Further, through spe- 
ciflc heat measurements on various samples of Yb2Ti207, 
we reveal evidence for sample dependence in both single- 
crystalline and powder samples, which we posit affects all 
samples reported on in the literature to date. Through 
neutron powder diffraction and Rietveld analysis, we 
show that the sample-dependence is not due to an ob- 
vious structural distortion, but rather to small levels 
of random disorder or strain within an average ordered 
structure. 



II. EXPERIMENTAL METHOD 



Ti02 in stoichiometric ratio and annealing at 1200°C for 
24 hours, with a warming and cooling rate of 100°C/h. 
The speciflc heat of this polycrystalline material was 
measured for comparison to the single crystals. 

Single crystal time-of-flight neutron scattering mea- 
surements were carried out using the Disk Chopper Spec- 
trometer at the National Institute of Standards and 
Technology (NIST), with 5 A incident neutrons, giving 
an energy resolution of 0.09meV. Triple-axis neutron 
scattering experiments were also performed at NIST, 
using SPINS. For the triple- axis measurements, a py- 
rolytic graphite (PG) monochromator provided an inci- 
dent energy of 5meV, and an energy resolution of ~ 8E 
= 0.25meV. Elastically scattered neutrons were selected 
using flve flat PG analyzer blades. Due to the aforemen- 
tioned energy resolution, the nominal elastic scattering 
in fact integrates over much of the quasi-elastic compo- 
nent of the diffuse scattering found in Ref. [l3|. A cooled 
Be fllter was used to remove higher order wavelengths 
from the scattered beam, and was followed by an 80' 
radial collimator. 

A neutron powder diffraction experiment was per- 
formed on two samples of Yb2Ti207, using the NDPF 
instrument at the Lujan Neutron Scattering Center. One 
sample studied was a polycrystalline rod (i.e. "sintered 
powder"), and the other was a single crystal grown in an 
identical manner to the sample used in the single crystal 
neutron experiments. Both samples were subsequently 
crushed using a Fritsch Pulverisette 2 mortar grinder, 
for 20 minutes per sample. 

The heat capacity measurements were performed us- 
ing the quasi-adiabatic method (see Ref. |30| for de- 
tails) with a 1 kr^ Ru02 thermometer and 10 Ml heater 
mounted directly on the thermally isolated sample. A 
weak thermal link to the mixing chamber of a dilution re- 
frigerator was made using Pt-W (92% Pt, 8% W) wire for 
the smaller single crystal piece, A (142. 3mg) and poly- 
crystalline sample (24.37 mg), and yellow brass foil for 
the 7.0515 g crystal, B, used in the neutron scattering ex- 
periments. The time constant of relaxation provided by 
the weak link was several hours, much longer than the 
internal relaxation time of the samples, minimizing ther- 
mal gradients and ensuring that the sample cooled slowly 
into an equilibrium state. The addenda contributed less 
than 0.1% to the speciflc heat of the system. 



Single crystalline Yb2Ti207 was prepared at McMas- 
ter University using a two-mirror floating zone image 
furnace. It was grown in 4 atm of oxygen, at a rate of 
5 mm/h, the growth procedure resembling Ref. [29|. The 
single crystal that was the subject of these neutron mea- 
surements was the same one studied previously ^^. The 
speciflc heat measurements were performed on this and 
one other crystal prepared under identical conditions. 
The starting material for these growths, a pressed poly- 
crystalline sample, was prepared by mixing Yb2 03 and 



III. RESULTS 

A. Diffuse Magnetic Neutron Scattering 

Figure [1] (a) shows the temperature dependence of the 
intensity of elastically scattered neutrons at two Q points 
using SPINS. One point is near the fee zone center and 
structurally allowed (111) Bragg peak (Q1.07 = [1-07, 
1.07, 1.07]), while the other is near the fee Brillouin zone 
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FIG. 3: a) Full width at half maxima (FWHM), and b) inte- 
grated area extracted from the Lorentzian fits to the trans- 
verse scans shown in Fig. 2a). c) The range of ^t within the 
kagome layer. Above 400mK (black) the correlations are re- 
stricted to a single kagome star unit. Below 265mK (blue), 
the transverse correlations grow only to encompass the next 
near neighbors in the kagome layer, d) The pyrochlore lattice 
decomposed into kagome and interleaved triangular layers. 
Below 400mK, the longitudinal correlations extend only be- 
tween a single KTK unit. 



boundary (Q1.54 = [1.54, 1.54, 1.54]). Figure [J (b) lo- 
cates these two points relative to the rod of diffuse scat- 
tering along [111] and the Brillouin zone boundaries, us- 
ing reciprocal space maps of the elastic scattering at 30 
mK and 500 mK measured with DCS. These two inten- 
sities track each other at relatively high temperatures, 
but they separate below ~ 400 mK, with the scattering 
near the zone centre (at Q1.07) growing more strongly 
before leveling off at the lowest temperatures. The inset 
of Figure [TJa) shows the difference between the elastic 
intensities at Q1.07 and at Q1.54, and this resembles con- 
ventional order parameter behavior with a phase transi- 
tion near 400 mK. 

A longitudinal scan of the elastic scattering along the 
rod (i.e. along the [H,H,B[] direction) shows that the low 
temperature buildup of intensity occurs near the (111), 
but not (222), zone centre (Figure 2a). Below 400mK, 
the scattering peaks up near (111) and the the Q de- 
pendence of the longitudinal scattering can be fit to an 
Ornstein-Zernike (Lorentzian) form, from which an in- 
terplane, or 'longitudinal', correlation length ^l can be 
determined. We find ^l = 3.2A at 60mK, implying cor- 
related spins over ~ 6. 4 A in the [111] direction. This 
indicates that even at the lowest temperatures, 3D corre- 
lations do not extend beyond a KTK stack, whose spatial 
extent is 5. 8 A. Above 400mK, this interplane correlation 
vanishes, leaving the Yb^+ form factor^^ as the only con- 
tribution to the Q dependence of the diffuse scattering 



along [111] as expected for truly 2D spin correlations. 

The scattering normal to the rod probes spin corre- 
lations within the kagome and triangular planes. Figure 
2(b) shows scans across the rods centered on Q1.07 and 
Q1.54. Representative scans at T=60 mK and T=540 mK 
are shown. These data sets were fit to Lorentzian line- 
shapes, allowing the transverse correlation length and 
integrated area of the scattering to be determined as a 
function of temperature. The resolution widths appro- 
priate to the spectrometer are shown as the horizontal 
bars in Fig. 2(b), and are negligibly small compared to 
the FWHM of the rod scattering. The results of the 
Lorentzian fits are plotted in Fig. 3 (a) and (b). The 
integrated areas at the two Q positions differ only for 
T<400mK, (Fig. 3 (b)) consistent with the behavior of 
the peak intensities along the rod shown in Fig. 1(a). 

Meanwhile, the FWHM (Fig. 3a), which is inversely 
related to the transverse correlation length ^t, does not 
depend on the longitudinal position along the rod at any 
temperature. Below 265mK, ^t is constant at 9.6A and 
it decreases to 7. 8 A above 400mK, leaving a correlated 
area corresponding roughly to one kagome 'star', which 
has a diameter of 15. 4 A. Combined with ^l^ above 400 
mK, we conclude that for T^ Qcw the spin system is 
characterized by uncorrelated kagome stars. Figure 3(c) 
summarizes the size of the spin correlated regions in the 
kagome planes both below 265mK and above 400mK. In 
the ground state, which stabilizes below 265mK, short 
3D spin correlations form along [111], encompassing a 
single KTK stack of width ~ 6A, as illustrated by the 
blue region in Fig. 3(d). 



B. Quasi-Elastic Neutron Scattering 

Inelastic magnetic scattering can inform on the dy- 
namics associated with this exotic disordered zero field 
ground state. Time-of-flight neutron scattering data is 
shown in Fig. |4] at both 30mK and 500mK, and at po- 
sitions both along the [111] rod of scattering and well 
removed from the rod of scattering. These positions are 
located in the reciprocal space map of the elastic scat- 
tering shown in the Fig. |4] inset. Consistent with earlier 
work, this shows the QE scattering along the [111] direc- 
tion extends to ~ 0.3 meV, but that the spectral weight 
is depleted at energies less than 0.2 meV, on entering 
the ground state. Hence, a dramatic slowing down of the 
spin fluctuations is observed, consistent with the ear- 
lier Mossbauer and /iSR results, ^^ as the 3D correlated 
ground state forms. Application of a [110] magnetic field 
in excess of ~ 0.5 T, induces a polarized, ordered state, 
and completely eliminates the diffuse scattering. Thus, 
the data in Fig. 3 (c) at 30 mK and 5 T serves as a 
measure of the background. 

Comparison of the of zero-field quasi-elastic diffuse 
scattering to the non-magnetic background revealed in a 
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FIG. 4: Energy dependence of the diffuse scattering above 
and below Tc, both on and off the rod of scattering. Thick 
hnes are guides to the eye. Inset: range of binning for the 
on- and off-rod positions. On-rod: HHH=[1.2,1.6], HH-2H= 
[-0.1,0.1]. Off-rod: HHH=[1.2,1.6], HH-2H=[0.3,0.5]. 
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magnetic field of 5T makes it clear that although diffuse 
scattering is organized into rods of scattering along [111], 
there is appreciable diffuse, inelastic scattering through- 
out the Brillouin zone. This "off-rod" scattering has a 
relatively flat energy dependence out to 0.6meV. Fur- 
thermore, unlike the rod diffuse scattering along [111], 
the "off-rod" diffuse scattering shows no temperature 
dependence between 500mK and 30mK. This implies 
that a dynamic and essentially uncorrelated component 
of the spins persists in the zero-field ground state. A 
similar situation is encountered in the classical spin ice 
Ho2Ti207, which retains appreciable elastic diffuse scat- 
tering throughout the Brillouin zone at very low temper- 
atures, indicating static zero-field spin correlations that 
encapsulate at most a single tetrahedron^ 



C. Specific Heat 

The specific heat provides insight into the extent to 
which the entropy changes at temperatures relevant to 
the diffuse scattering. Figure [5] shows specific heat re- 
sults of poly- and single-crystalline samples, including 
the sample from the neutron scattering study. There is 
significant sample dependence of the low temperature 
Cp anomaly, with at least one feature occurring between 
150mK and 265mK in all samples. The sharpest anomaly 
is observed in our powder sample, at a temperature that 
seems to be the upper limit for all others, T=265mK. 
Significantly, our powder Cp peak is sharper by an order 
of magnitude, and occurs at a higher temperature than 
both the original powder Cp data by Blote et al^ as well 



FIG. 5: Examples of low temperature specific heat in 
Yb2Ti207. Significant sample dependence is observed in 
both powders and single-crystals. The powder sample pre- 
pared at McMaster University shows the highest temperature 
(265mK) and sharpest anomaly. The neutron scattering sam- 
ple, B, shows a sharp peak at 265mK similar to the powder, 
and a broad, low temperature feature similar to crystal A. 
None of the samples studied had any feature at 400mK. 



as the powder Cp data by Dalmas de Reotier et al, which 
displays an anomaly at 250mK with a peak height of ~9 
J/K mol Yb^ The neutron scattering sample B exhibits 
a sharp peak at 265mK, but also has a broad feature as 
seen in sample A. 

The specific heat does not show any indication of a 
feature at 400mK. This indicates that the buildup of 3D 
correlations, as well as the small increase in transverse 
correlation length, produces, at most, subtle changes 
in the entropy. This could be another manifestation of 
the persistently short correlation lengths at all temper- 
atures. However, an important question remains - what 
is the physical significance of the Cp features at 265mK 
and 200mK observed in the neutron scattering crystal? 
Drawing on the large change in spin fluctuation fre- 
quency observed by Hodges et al at 240mK in their pow- 
der sample^ we should expect a change in dynamics at 
these temperatures. As discussed in relation to Fig. 4, 
we do observe a change in dynamics between 500mK and 
30mK but the intermediate temperature regime awaits 



further study. 



D. Structural Information 

To investigate structural differences between the sin- 
gle crystals as grown by the method described in Sec- 
tion II, and the polycrystalline material which was used 
as the starting material for these growths (the "sintered 
powder"), we turn to neutron powder diffraction results. 
After a successful growth, a single crystal was crushed 
back into polycrystalline form in order to facilitate di- 
rect comparison to the sintered powder. Figure [6] shows 
the two neutron diffraction patterns from the powder 
and crushed single crystal at 250K. The data is nor- 
malized to the incident time-of-flight spectrum, and cor- 
rected for detector efficiency using a vanadium normal- 
ization. A Rietveld refinement was performed on both 
patterns using a fully ordered model of the pyrochlore 
structure (see, for example, Ref. 32 for details of this 
structure). The results indicate that both the powder 
and crushed crystal are described extremely well by the 
fully ordered model, with Rp values of 2.16% and 2.83% 
respectively. Some differences can be observed between 
the two samples, namely that the crushed crystal has a 
larger lattice parameter {axtai = 10.02006(2) compared 
to apow = 10.01322(1)) as well as greater widths of the 
Bragg peaks, indicating lattice strain. Site substitution 
between Yb^+ and Ti^+ was found not to be important 
for either sample (refinement gives less than 1% substitu- 
tion in both cases). A detailed, temperature-dependent 
Rietveld refinement as well as Pair Distribution Function 
analysis is forthcoming, but these preliminary results al- 
ready indicate that there are only subtle structural dif- 
ferences between the two forms of Yb2Ti2 07, pointing 
to random disorder or lattice strain as the root cause for 
the differences in specific heat. 



IV. CONCLUSIONS 

In summary, we have explored, using diffuse and in- 
elastic neutron scattering, the transformation of spin cor- 
relations as Yb2Ti207 enters its exotic ground state at 
zero field. Within this ground state, we observe slow, 
short-range 3D correlated fluctuations that extend over 
single KTK stacks. In contrast, at temperatures charac- 
teristic of Qcw and above, the system is comprised of de- 
coupled kagome stars, with 2D spin correlations confined 
to single kagome units. The temperature T* :^ 400mK at 



which correlations lose their 3D nature and also shrink 
in the 2D kagome plane is found to be much higher than 
any of the temperatures Tc < 265mK marking anoma- 
lies in Cp. Random structural disorder or lattice strain 
may account for the variability in Tc and the sharpness 
of the associated anomalies, but would seem unlikely to 
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FIG. 6: Neutron diffraction data from two samples of 
Yb2Ti207 taken at 250K: a sintered powder (top) and a 
crushed single crystal (bottom) . The result of Rietveld analy- 
sis, using a fully ordered model for Yb2Ti2 07, is shown in red. 
The difference between the measured and calculated profiles 
is in black. 



be the factor capping the extremely short range of cor- 
relations in the ground state. Both T* c:± 400mK and 
the highest Tc = 265mK observed are an order of mag- 
nitude smaller than the predicted mean-field transition 
temperature, 3.2Kj^ It is clear that the transition to 
the exotic ground state, whether this be characterized 
by Tc and/or T*, is strongly suppressed by geometrical 
frustration, quantum fluctuations, or both. 
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